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On examining the stability of a Wigner crystal in an ionic dielectric, two competitive 
effects due to the polaron formation are found to be important: (i) the screening of the 
Coulomb force X/e s r, which destabilizes the crystal, compensated by (ii) the increase 
of the carrier mass (polaron mass). The competition between the two effects is care- 
fully studied, and the quantum melting of the polaronic Wigner crystal is examined by 
varying the density at zero temperature. By calculating the quantum fluctuations of 
both the electron and the polarization, we show that there is a competition between the 
dissociation of the polarons at the insulator-to- metal transition (IMT), and a melting 
towards a polaron liquid. We find that at strong coupling, a liquid state of dielectric 
polarons cannot exist, and the IMT is driven by the polaron dissociation. Next, taking 
into account the dipolar interactions between localized carriers, we show that these are 
responsible for an instability of the transverse vibrational modes of the polaronic Wigner 
crystal as the density increases. This provides a new mechanism for the IMT in doped 
dielectrics, which yields interesting dielectric properties below and beyond the transition. 
An optical signature of such a mechanism for the IMT is provided. 



1. Introduction 



In a recent letter to Nature, Phillips et aL0 pointed out that a superconducting 
state could be expected close to the instability of a 2D Wigner Crystal (WC) of 
electrons. Previously, Takadaa had mentioned such a possibility for a 3D WC. In 
both cases, the pairing mechanism is due to correlation between electrons, mediated 
by the long range Coulomb interaction. More explicitely, pairing occurs through 
the Coulomb holes which form in the vicinity of the melting transition. Translated 
into the language of dielectric constants, as was first recognized by Bagchfl for the 
3D WC, it means that the static dielectric constant is negative. In his original 
paper, Bagchi also extended his results to an insulating lattice of Lorentz dipoles. 
He found that, on increasing the dipole density, the Lorentz lattice undergoes a 
phonon instability due to the dipolar interactions. At the same time, the static 
dielectric constant becomes identical to the dielectric constant of an ordinary WC, 
and thus can be negative. In other words, an ordinary WC appears to be nothing 
but a Lorentz lattice of dipoles just on the verge of instability. This could provide 
a natural pairing mechanism for (hypothetical) free electrons moving through such 
a Lorentz lattice close to its critical point. Quite generally, the question of the 
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sign of the static dielectric constants has been addressed by many authors. Cohen 
and AndersonQ argued on physical grounds that such a negativejign would be 
impossible. Later, the "rassian school" , in particular Dolgov et aln, following the 
previous work of MartinQ, pointed out that a negative sign of the dielectric constant 
would not contradict any stability criterion. The only required condition is that 
e (0, 0) be strictly positive: a system can have e (k, 0) < for any k > 0, and still 
be stable. 

We have developed!] such a physical picture in order to describe the Insulator- 
to-Metal transition in ionic dielectrics, taking into account the formation of large 
polarons when introducing the doping charges. These polarons, which are ordered 
in a crystallized state at low densities, behave more or less as Lorentz dipoles, 
and can yield some of the interesting physical properties already mentioned by 
Bagchi. In order to understand our point of view, it is useful to connect our studies 
to the recent "fight" in Physical Review Letters of Chakraverty et at against 
the (bi)polaronic superconductivity scenaridlj. These authors did not find the 
bipolaronic scenario suitable for the interpretation of the experimental data in high- 
Tc, a point of view which seems to be fully shared by Andersonu. In fact, we 
would like to underline that the basic building block of the (bi)polaronic scenario 
is the existence of a quantum liquid of (bi)polarons. As we will explain later, 
this turns out to be impossible at strong electron-phonon coupling, owing to the 
drastic role of long-range Coulomb forces in the many-polaron problem. Indeed, 
treating self-consistently a finite density of interacting polarons is a much more 
complicated problem than considering only one or two polarons. Extensions to finite 
densities usually neglect the Coulomb intera ctio n between (bi)polarons, as in the 
framework of the random phase approximations. As we will see, this approximation 
is quite dubious at strong polar electron-phonon coupling, while it should remain 
reasonable at small and intermediate couplings. Long-range Coulomb forces are 
responsible for both the polaron formation (if one only considers the coupling of 
the electrons to the polar mode of the crystal, as it is the case in all our studies), 
and for their mutual interactions. In such systems, electron-phonon and electron- 
electron interactions cannot be decoupled and must be treated on the same ground, 
otherwise the collective behavior of polarons at finite density can be completely 
missed. To tackle the many-polarom_riroblem, we have considered a system at low 
density. Starting from the remarkQa that two polarons at a distance d_cfipel as 
l/e s d if the dielectric constants satisfy £oo/e s >~ 0.1 (as it is well-knowr£3, there 
is no bipolaron formation in that case), at low densities, the crystallized state is 
necessarily the ground-state. 

Let us consider a WC, consisting of a system of electrons plus a jellium (assumed 
to be rigid) of positive charge with the same density. We can imagine that this WC 
is inserted into a host polar material, characterized by the two limiting dielectric 
constants e s , and by the frequency of the longitudinal optical mode ojlo- If 
one turns on the coupling between the WC and the polar phonon mode of the host 
material, two competitive effects arise. The first one corresponds to the screening 
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of the Coulomb interaction between the electrons, which is reduced from l/e^r to 
l/e s r. This screening is due to the host lattice polarization which surrounds the 
electrons, which is nothing but the polaron cloud. The other important effect is 
the increase of the carrier mass: when the coupling with the lattice is turned on, 
the original band mass m* becomes Mp, the polaron mass, which can be large, 
especially in the strong coupling limit for which the polarons are well-formed. Now, 
the crystallization is due to the overcoming of the kinetic energy l/M p R 2 s by the 
Coulomb repulsion l/e s R s , where R s is related to the density n as usual: n = 
We see from this simple argument that the polaronic crystallized state 
can be stable up to fairly high densities in spite of the screening of the Coulomb 
forces. Furthermore, as given by standard polaron theory, the mass Mp becomes 
infinite as the phonon frequency tends to zero. In that case, i.e. when ojlo ~ * 0, 
the kinetic energy is always negligible with respect to the Coulomb repulsion, and 
melting towards a liquid state of polarons is clearly not possible. 

In our studies, we have focused on the quantum melting (at zero temperature) 
of such a polaronic Wigner crystal (PC). We have found that the two competi- 
tive effects above give rise to two different scenarios for the IMT, which can be 
easily understood. When the electron-phonon coupling vanishes, we are left with 
an ordinary WC of electrons. As the density is raised, the quantum fluctuations 
of the localized particles increase, and when they reach a certain magnitude, the 
crystal melts towards a liquid state of electrons. The effect of the host lattice is 
negligible in that case. In the opposite limit lulo ~^ 0, the electron-phonon corjb 
pling becomes infinite. Its strength is given by the usual dimensionless parametert-J 

1/2 

a = (m* /2h 3 LULo) e 2 /e, where 1/e = 1/Eoq — is the effective dielectric 
constant responsible for the polaron formation. In this limit, each electron is self- 
trapped in a polarization potential-well which is frozen. As it is well-known, the 
latter has a coulombic nature and behaves as 1/er at large distance. Since the po- 
larization cannot move, the polarons, i.e. the electron plus the polarization, cannot 
go towards a liquid state. The IMT is thus driven by the dissocation of the polarons 
and the screening of the polarization potentials by the liberated electrons. We are 
in a situation which closely resembles the usual Mott transition. On the basis of 
the same argument as MottB, the critical density can be estimated to be: 

nj/ 3 (e/eoo )i? P « 0.25, (1) 

where Rp is the bound-state radius, which is nothing but the polaron radius given 
at strong coupling by Rp ~ 3.2e(m* /m e )ao- Here, ao is the usual Bohr radius 
(0.53A), and Rp corresponds to a hydrogenic wave-function. 

Starting from the limiting cases presented above, it is now important to ex- 
tend the results to any finite electron-phonon coupling (i.e. to any finite phonon 
frequency), and examine which mechanism is responsible for the melting of the 
polaron crystal. At this stage of the theory, we have developed two steps. 
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Figure 1: (a) The fluctuation of the polaron (5R 2 ) 1 / 2 / R s (bold line) and of its 
internal degree of freedom (dr 2 ) 1 ' 2 / R s (dotted line) versus n at a — 10. (b) The 
polaron fluctuation (SR 2 ) 1 ^ 2 / R s versus n for different a. Physical parameters are 
£co = 5, e s = 30 and m* = 2m e . The density parameter is simply n — I .6 • 10 24 /rf . 

2. Quantum melting: polaron breaking versus polaron liquid 

First, we have modeled the situation in the PC by considering a single polaron 
localized in a positively charged jellium sphere with density n — 3/AitR^. In this 
approach, which is a sort of mean-field many-body calculation (see for example 
the first papers by Wigner on the electron crystallization) , dipolar interactions and 
exchange effects are neglected. However, one obtains reliable results on the crys- 
tal energy and ground state properties, provided the_ density is low enough (large 
R s ). As it was done by Nozieres et al. for the WCEJ3, we have applied the Linde- 
mann criterion to this model. These calculations allow us to understand, at least 
qualitatively, the route that the crystal will choose to melt: polaron dissociation 
versus polaron liquid state. Of course, this method just gives access to the melting 
scenario, not to the ground-state beyond the transition. 

The application of this criterion to the PC is more complicated than for an 
ordinary WC, owing to the composite nature of the polarons. To be able to delocalize 
the polarons towards a liquid state, the quantum fluctuations for each polaron 
moving as a whole, i.e. the motion of each electron together with its surrounding 
polarization, must become large in comparison with R s . On the contrary, if the 
relative fluctuations of the electron with respect to the polarization become large 
the polarons break apart (dissociation) . The Feynman treatment for the polaronEZI, 
which relies on path integral calculations, yields a natural way to evaluate such 
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quantities: since the polarization field in Feynman's approach is replaced by a rigid 
particle with coordinate X, one easily calculates both the quantum fluctuations of 
the center of mass R — (m* x+ M X) / \m* + M) and the relative coordinate r — x~X 
(x is the electron coordinate). Thus, we now have two different Lindemann criteria, 
which correspond to two different melting scenarios: 

1/2 

(i) (SR 2 ) /R s > 1/4 which allows the melting of the crystal towards a liquid 
state of polarons; 

1/2 

(ii) (i5r 2 ) / R s > 1/4 for the polaron dissociation. 

Both quantities are calculated by the same path integral variational proceduro30 

The frequencies of the two degrees of freedom r and R, respectively Uint and 
uJext, were calculated in reft together with the ratios (i) and (ii) (see Fig.[j]) for 
different a. From the basic observation that a polaron is a bound state of an 
electron plus a phonon cloud, it was pointed out that the frequency of vibration of 
the polaron as a whole is physically limited by the phonon frequency (uo ext < ojlo)- 
In other words, at strong electron-phonon coupling (or equivalently small ojlo)i the 
lattice polarization cannot dynamically follow the increase in kinetic energy induced 
by the doping, and the quantum fluctuations are transferred to the internal degree 
of freedom, breaking the polaron apart. For that reason, the polarons are expected 
to dissociate at the transition for sufficiently high a (a > a* ~ 7.5), rather than 
melting to a polaron liquid. 

We have also calculated the same quantities for an anisotropic polaron. By 
anisotropic, we mean that the band mass of the electron is different in the xy plane 
than in the z direction, by means of the phenomenological parameter 6 = m xy /m z . 
However, the polarization is assumed to preserve its three-dimensional character. 
This is the case Jbr cuprates, as we have already dicussedQ, on the basis of Chen et 
al. experimental^. When 6 = 1, we recover the usual 3D polaron, while for 6 = 0, 
we have some kind of 2D polaron (the electron motion is restrained to the xy plane). 
Table 1 summarizes the numerical results obtained for two different band masses 
m* = m e and m* = 2m e , for both 6=1 and 6 = 0. 





a* n c (m* — m e ) n c (m* = 2m e ) 


isotropic (6 = 1) 

(fully) anisotropic (6 = 0) 


7.5 6-10 ls 5-10 iy 
4.4 8 • 10 19 6 • 10 20 



Table I. Critical density n c (cm 3 ) and coupling a* (see text) 



The important qualitative fact that we can conclude from the numerical data 
is the following: what happens for a > 7 in 3D, will happen for lower couplings 
(a > 4) in the anisotropic case. In other words, the polaron dissociation may be a 
physical possibility in the cuprates for describing the IMT mechanism, at view of 
the fairly high value of the polar electron-phonon coupling (Chen et alr3 estimated 
a ~ 5 for La 2 Cu04). 
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3. Effect of the dipolar interactions: Bagchi's instability 

Let us now concentrate on the melting in the 3D case, and restrain ourselves 
to the limit a — > oo (lolo — > 0). In this regime, the lattice polarization around 
each electron is frozen, and we can easily include the dipolar interactions that were 
previously neglected. The electrons vibrate in a potential-well due to the static 
polarization, the jcllium sphere and the electrons on the other sites. Assuming 
that these are localized on the sites of a Bravais lattice, the crystal of polarons can 
be viewed as a Lorentz lattice of dipoles. Each individual dipole vibrates at the 
frequency uji nt above mentioned. BagchJl has already calculated (classically) the 
phonon dispersion and the dielectric properties of such a Lorentz lattice. It can be 
shown! that the vibrational spectrum (k, A) of such a lattice (k is the wave vector 
and A = 1, 2, 3 is the phonon branch index) satisfies the following relations: 



n 2 (k, A) = u? (k, A) + w 



pol 7 



(2) 



where we have defined ui„„ 



-w 2 /3e c 



47rne 2 /m* being the usual plasma 
nothing but the 



1 @ is 
IS: y 



EL w ( k . A ) 2 



and 



frequency (defined without any dielectric constant), uj (k, A) in (2 
phonon spectrum of a usual WC. It satifies the Kohn sum rule 
ojp/soc- At long wavelengths, there is one optical mode (u (0, opt.) = oj p / \fe~ a 
two acoustical modes. From these properties of ui (k, A) and using (2), it comes out 
immediately that the long wavelength transverse modes of the polaron crystal become 
unstable when us po i < 0, i.e. when u>i nt < lu p /\/3£oo- The critical density for which 
w po ;(n c ) = 0, based upon the calculation of u>i n t in section 2, is n~ ag ' as 5 • 10 20 cto~ 3 
for m* = 2m ei e s = 30 and Eqo = 5, which is slightly higher than the value obtained 
through the Lindemann criterion. 

The (longitudinal) dielectric constant e (k, oj) of the Lorentz lattice has also been 
calculated by Bagchi. 



1 



1 



e(k,w) e c 



w 2 /e 



50 E 



( k • e k,A) 



t [^(k, A) + ^-a,2 



(3) 



where ek,A is the polarization wavevector of the phonon mode k, A. It can be easily 
evaluated both in the static and in the long wavelength limit. In the latter case, we 
have: 

Up/ £qo 



e(0,w) 



1 - 



^ 2 ~ ^pol 



(4) 



which is sketched in Fig. 2. a. The pole in the dielectric function, located at lu po i 
(i.e. the polaron peak in the optical conductivity), is shifted to lower energies 
as the density increases, and tends to zero on approaching the IMT, as shown in 
Fig.2.b. Strictly speaking, since ojlo is never zero, but assumed to be small, to po i 
should rather saturate at higher density. The corresponding behavior of the optical 
conductivity ct(cj) = lm[e(0, cj)]cj/47t and energy loss function Im[— l/e(0, u)] at 
T = is sketched in figure 3. The polaronic peak shift seems to have been recently 
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1 2 3 4 5 6 
n(10 20 cm- 3 ) 

Figure 2: (a) sketch diagram of the long wavelength dielectric function, (b) shift of 
the polaron peak and crystal instability in the static limit (ojlo = 0), for m* = 2m e . 



observedcJ in the cuprates. Even though our calculations are strictly valid only in 
the limit ujlo ~ * 0, one can however argue that for a > a* (above mentioned), the 
results - in particular the existence of Bagchi's instability - remain qualitatively cor- 
rect, since in that case the polarization cannot easily follow the electron vibrations. 
New kinds of collective modesEil corresponding to the mixture of polarization waves 
with electron vibrations can also be expected, as was pointed out by Bozovi<x3, 
which should have an optical signature. 

As mentioned in the introduction, the static dielectric constant also exhibits in- 
teresting properties in this limit. From (3), it can be seen that at the critical density 
(ujpoi — 0), owing to the sum rule, the phonon dispersion obeys to (k, A) < oj p / ^j£oa- 
Therefore, we necessarily have e (k, 0) < for any k > 0. This indicates the pos- 
sibility of having a net attraction between free carriers coexisting with localized 
polarons beyond the transition. This possibility opens up very interesting ques- 
tions for both a superconducting scenario, and unusual properties of the normal 
phase beyond the IMT. Pseudo-separation between kxcalized (heavy) polarons and 
free (light) electrons has already been mentionecu'OSo both theoretically and 
experimentally. We believe that Bagchi's instability can be one of the keys of this 
pseudo-separation. Evidently, at the critical density, some polarons must be de- 
stroyed. However, it is not necessary to destroy all the polaron lattice, and we just 
need to fulfill the stability condition on the vibrational modes il 2 (k, trans.) > 
(i.e. ripoiarons < nf a9 -). On the other hand, there is some energetical argument 
in favor of a mixed phase of localized polarons with free electrons. The polaron 
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(o/to pol (0) <a/GV,,(0) 



Figure 3: (a) optical conductivity and (b) EELS spectrum, as calculated by in- 
troducing an arbitrary Lorentzian broadening in eq.(4). Labels 1-5 stand for the 
density in units of 10 20 cto~ 3 . Physical parameters are the same as in Fig.l, and 
n c « 5 • 10 20 cm~ 3 . 

formation energy being negative, it costs energy to destroy them. A compromise 
might be achieved (which remains to be proved) by the coexistence of a polaron 
lattice with density nf a9 \ and quasi-free electrons, provided the mean distance be- 
tween free electrons (corresponding roughly to the screening length) is larger than 
the distance between localized polarons. However, above a certain density (say 
rtfree ~ n^ ag ), all the polarons must disappear due to a global screening effect. 
This last situation could correspond to the overdoped region in the cuprates, while 
the coexistence between localized polaron and free electrons, to the underdoped 
region. 

4. Conclusion 

The studies proposed in this paper (and previous references) are currently being 
extended in several directions, essentially at iinite temperature and finite phonon 
frequency, including the dipolar interactionsc!! The possibility of a mixed state 
of crystallized polarons coexisting with quasi-free carriers is important and must 
be clarified. Another problem concerns the possibility of having a crystal of bipo- 
larons instead of polarons. This is based on the recent results of Moulopoulos and 
Ashcroft£3 for the ordinary WC, where they have proved that at intermediate den- 
sities, a paired electron crystal is more stable than the single electron crystal. In the 
polaronic case, this is still an open question. Nevertheless, it would not change the 
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basic result of our studies: at strong coupling, the IMT is driven by the dissociation 
of some localized (bi)polarons and cannot yield a liquid state of (bi)polarons. 

In conclusion, let us underline that what we have called the Bagchi instability, 
can also be supported by the existence of localized dipoles which have nothing to do 
with polarons. This remark comes in mind if one considers another family of doped 
oxides, the bismuthates (K x Ba\- x BiOz), where it has been recently suggested that 
dipolar interactions could be of importanceEj. 

Part of this work has been presented at the International Conference STRIPES 
in Rome (2-4 june 1998), and took its origin from very long discussions that one of 
us (P.Q.) had with A. Bianconi during the year 1993, on the possible crystallization 
of polarons in the cuprates. This work received financial support from the European 
Commission (contract no. ERBFMBICT 961230). 
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